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Anomaly-induced charges in baryons 
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We show that quantum chiral anomaly of QCD in magnetic backgrounds induces a novel structure 
of electric charge inside baryons. To illustrate the anomaly effect, we employ the Skyrme model for 
baryons, with the anomaly-induced gauged Wess-Zumino term ~ (tyo + multi-pion)i?- B. Due to this 
term, the Skyrmions giving a local pion condensation ((tto + multi-pion)) / necessarily become a 
local charge source, in the background magnetic field B ^ 0. We present detailed evaluation of the 
anomaly effects, and calculate the total induced charge, for various baryons in the magnetic field. 



I. INTRODUCTION 

The chiral anomaly is one of the central concepts in 
QCD, and it manifests nature of quantum field theories 
in an explicit way in our hadronic world. As the chiral 
anomaly is essentially coupled to electromagnetic sector 
since the electromagnetism is a part of the chiral symme- 
try, the introduction of nontrivial electromagnetic back- 
grounds should add a good flavor of physics onto the 
chiral anomaly. In this paper, we report an interesting 
new effect induced by the chiral anomaly, for baryons in 
a background magnetic field. 

Our finding is that baryons in a constant magnetic 
background acquire additional electric charge distribu- 
tion due to the chiral anomaly. The result that this 
would generate even a total net charge is quite surpris- 
ing, but the mechanism is quite simple. It is well-known 
that Wess-Zumino- Witten (WZW) term @, Q actually 
captures the chiral anomaly in terms of the hadronic de- 
grees of freedom. In particular, this term serves as a 
manifestation of the famous tto — >> 2j decay. Now, any 
baryon carries a cloud of pions around it, and so it is 
a source of the pions. Once we replace one of the two 
7's in the Wess-Zumino- Witten term by the background 
magnetic field, we immediately see that the baryon can 
be a source of the electromagnetism (another 7), i.e. the 
baryon can have an additional charge structure due to 
the chiral anomaly and the pion cloud. The schematic 
picture of this mechanism is illustrated in Fig. [TJ 

In this paper, we explicitly demonstrate this mecha- 
nism in detail, with a help of a concrete model of the pion- 
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FIG. 1: A schematic figure for electric charge generation of a 
nucleon. In electromagnetic backgrounds, i.e., F^ v / 0, the 
chiral anomaly generates an additional coupling to the gauge 
fields A u . 



cloud picture of the baryons, the Skyrme model [1] . In 
the Skyrme model, baryons are given as a solitonic object 
made of a local pion condensate (tt(x)) 7^ 0. Plugging the 
Skyrme solution to the Wess-Zumino- Witten term, it can 
be shown that the magnetic field background can induce 
a novel charge structure inside the baryon (Skyrmion). 

In particular, we give an argument that the total 
charge can also be generated, and resultantly the Gell- 
Mann-Nishijima formula for baryon charges can be cor- 
rected under the magnetic field due to the anomaly, 



= e[h 



Nb 
2 



(1) 



Here in the modified formula, Q e is the electric charge of 
the baryon, 1% is the third component of isospin, Nb is 
the baryon number, and the new term Qanm is the charge 
generated by the anomaly and the background magnetic 
field. 

One may be suspicious on this generation of electric 
charges. However, for example in the renowned Wit- 
ten effect 0, [Bj] , monopoles are accompanied with electric 
charges, in the presence of the term. We may regard 
our WZW term as an analogue of the term for the 
Witten effect. In addition, the chiral magnetic effect 
0] in heavy ion collisions shares the same property too. 
So, it is fare to say that the generation of the electric 
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charge is not a unique feature of our investigation, but is 
a common feature among parity- violating effects. 

Quantum anomaly is literally quantum-mechanical, 
and thus is a tiny effect. However, when the coupled 
magnetic field is strong, this effect may be enhanced. So 
our physical motivation for this work is primarily ori- 
ented to the situation in which strong magnetic field is 
present with a finite density of baryonic matter. For this, 
one can come up with two important physical cases: one 
is a neutron star, at which neutrons are very dense and 
with a strong magnetic field, and the other is a heavy 
ion collision at which nuclei are smashed and a strong 
electromagnetic field is expected to be created instantly. 
In this paper, we do not go into these concrete cases. We 
concentrate on providing a basis for that, and in partic- 
ular evaluate in detail the anomaly WZW term with the 
quantized Skyrmions, under a constant magnetic field. 

The organization of this paper is as follows. In section 
III1 we provide a review of the Skyrme model and the 
WZW term, with a brief introduction to the Skyrmion 
solution. In section Hill we shall see explicitly that the 
background magnetic field generates an additional charge 
structure in the Skyrmions (baryons). We quantize the 
Skyrmion and evaluate the anomaly-induced electric cur- 
rent for an arbitrary baryon state. In section HV1 we eval- 
uate the multipole moments of the anomaly-induced elec- 
tric current and found a quadrupole, with a pion-mass de- 
pendence. In section [VJ we discuss possible other effect 
due to the background magnetic field on the baryon. In 
section [VlJ we evaluate classically the anomaly-induced 
charge for higher-charge (=multiple) Skyrmions. The fi- 
nal section is for our conclusion and discussions. Ap- 
pendix [A] is a study of the generated charge in a non- 
constant magnetic field. In Appendix [Bj we show that 
the induced charge is due to a multi-pion effect (i.e. a 
pion cloud) , and we compare our result with a point par- 
ticle description of baryons. The letter version of this 
paper is [10| . 



II. THE SKYRMIONS AND THE 
ANOMALY-INDUCED CHARGES 



As briefly described in the introduction, it is indeed 
almost straightforward to calculate the effect of the 
anomaly term for baryons in the presence of the mag- 
netic field background, once we adopt a concrete model 
of the pion cloud. Here, we first review the Skyrme model 
which realizes baryons as a condensation of the pions, and 
also review the gauged WZW term which manifests the 
chiral anomaly in QCD. 



A. The model 

1. The Skyrme model 

The chiral symmetry SU(N) L x SU(N) R acts on left- 
handed and right-handed quarks as 

<7L U L q L , q R U R q R: with U LiR e SU(N) L , R .(2) 

When the chiral condensate q R qh develops a non-zero 
vacuum expectation value by some non-perturbative ef- 
fects 

(qRqL) = -v s l N , with v = O(Aqcd), (3) 

the axial-part of the chiral symmetry is spontaneously 
broken as 

SU(N) L x SU(N) R -+ SU(N) L + R . (4) 

This gives rise to Nambu-Goldstone (NG) bosons, 
namely the pions, which takes value in the coset space 

SU(N) L xSU(N) R 
SU(N) L+R ' 

U{x) = exp (^^T^J , (a = 1, 2, • • • , N 2 - 1). (5) 

Here F n = 108 [Me V] is the pion decay constant and T a is 
a generator of SU (N) and we use the following standard 
normalization 

Tr[T a T b ] = h ab . (6) 
The chiral symmetry acts on the NG modes as 

U U L UUl (7) 

For later convenience, let us define left- and right- 
invariant Maurer-Cartan one-forms by 

L^U^d^U, R^ = d^UUl (8) 

These take their values in the algebra of SU(N) R and 
SU (N)l, respectively. The chiral symmetry acts on them 
as 

i/i -+ U R L^Ul -+ U L R^Ul (9) 

We can think of U as an effective low-energy field. Its 
effective Lagrangian of the leading order to 0(d 2 ) can be 
uniquely determined as 

F 2 

£ (2) = ^^[d^Ud^ + M%(U + -2)] 
F 2 

= -JLTi[-R^ + Ml(U + U^ -2)]. (10) 

Here M n stands for the pion mass M n = 137[MeV] and 
our metric is r]^ v = diag(+l, —1,-1,-1). By expanding 
L M and with respect to l/i^, one gets 

L M = 4i%^T a + 8ie abc ^f^T c + ... , (11) 
Rfi=4i ^T a -8ie abc ^ LTC + --- ■ (12) 
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Plugging this into C^ 2 \ one obtain a standard kinetic 
term of the pions and corrections, 

1 M 2 
£( 2 ) = -9/9V--iiY 

2 p 2 



9 M 2 



(13) 



We are interested in a topological soliton made by the 
pions in this work. The topological winding number is 
given by 



7T 3 (SU(N)) =N B eZ. 



(14) 



As will be shown, Nb is identified with the baryon num- 
ber via the WZW term. However, it is easy from a simple 
scaling argument that no topological solitons can survive 
from collapsing in the theory with C^ 2 \ So one needs 
higher derivative corrections to C^ 2 \ Therefore, we take 
a term of order 0(d A ) which is so-called the Skyrme term 



£ (4) 



1 



32e? 



Tv{[R^R„}[R»,R»}), 



(15) 



where we have used 

UL^U* = f/([/t^[/)[/t = R v m 



(22) 



The equation of motion of the Skyrme model is identical 
to the current conservation law if the pion mass is zero 

8^ = 0, or 8^ = 0. (23) 

When the pion mass is non-zero, the equation of motion 
becomes 



« = -^|^Tr[C7-f/t]. 



(24) 



The vector and axial conserved currents are defined by 

^ = i£+i| ft = &^&. (25) 



The vector SU (N) l+r is nothing but the isospin, so we 
write its conserved charge as 

I a = J d 3 x f v a = J d 3 x Tr[(f L + f R )T a }. (26) 



with e s being a dimensionless coupling constant. We will 
choose the parameter e s = 4.84 by following Ref. [l2j . 
Now we are ready to write down the Skyrme model with 
the right- invariant one form as 



£ = -^Tr [-R^ + M^(U + U^ -2)] 
' :Tv([R^R u }[R^R"}). 



32e- 



(16) 



A Noether current of SU (N)l can be obtained by per- 
forming a local and infinitesimal SU(N)l rotation 

SR^id^L. (17) 

Variation of the Skyrme Lagrangian is given by 

SU = Trjg (-FlR» + ^[R v , [R»,R»}^ 8^ .(18) 

Then the conserved current is given by 



ft = ±(tfR> t -±[ir,[R ll ,R v ] 

Similarly, the SU(N)r current takes the form 

1 



-F^ + ^[L V ,[L^,L V 



These currents are related by 



(19) 



(20) 



(21) 



2. Electromagnetic interaction 

Let us next take the electromagnetic interaction into 
account. For simplicity, hereafter, we concentrate on the 
minimal case with two flavors N = 2. Since the elec- 
tric charges of u and d quarks are 2/3 and —1/3 respec- 
tively, the NG modes are rotated under the electromag- 
netic U(l) as 



U ^ e ^eQ U ^eQ = e -ieT* Ue ieT* 



Q = il + T 3 ,(27) 
o 



where T 3 = r 3 /2 with the Pauli matrix r a . Thus the 
electromagnetic U(l) em is a subgroup of SU(2)l+r. 

Interactions of the NG modes and the electromagnetic 
fields are introduced by gauging the U (l) em C SU (2) l+r 
and replacing the partial derivative by a covariant 
derivative 

VpU = 8^U + ieA^[T 3 ,U}. (28) 

The left- and right- invariant one-forms are then replaced 

as 

R^R^= V^UU\ L^L^ = rfV^U. (29) 
Then the total Lagrangian can be read as 



1 

~4 
1 

32e 



F„ V F^ 



F, 



if 1 * 



-RpB? + (U + t/ f - 2) 



l ^{[R l _ k ,R v }[R lx ,R v \). 



(30) 



The classical equation of motion is derived by varia- 
tional method as before. One can easily obtain 



i F 2 M 2 



(31) 



One can also express the E.O.M. in terms of the right- 
invariant one- form L M by just replacing with L M . Note 
that since the electromagnetic charge Q breaks the chiral 
symmetry SU(2) L x SU(2) R /7(l) em , ]£ and 3r are 
not conserved currents. The Maxwell equation is given 

by 

= eft", ft v =Tr[T 3 (ji+f R )]. (33) 
Thus the electromagnetic charge is given by 



■jd 3 x~j 3 »=° 



el 3 



(34) 



3. WZW term and Chiral anomaly 

In order to describe the baryons in the Skyrme model 
whose fundamental degrees of freedom are mesons, we 
have to consider not only elemental particles but also 
topological excitations. Since the NG modes do not carry 
the U(1)b charges, we need to add an extra terms to 
the above Lagrangian. It is the so-called Wess-Zumino- 
Witten (WZW) term. With an electromagnetic field 
the WZW term in the N = 2 flavor model is given by 



[14] 



(35) 



with the baryonic current 



r B 



i 



; e^Tt[R u R p R a ] - --^e^dv (A p P a ) .(36) 



24tt 2 -l-^v— j g7r 2- 
Here we use e 0123 = — 1 and define 



P a = -Tr[ T 3 (L a + R a )}. 



(37) 



This baryon current is clearly conserved due to the anti- 
symmetric tensor e pupa . On the other hand, j£ appears 
to depend on gauge choices, at a glance. But this is not 
the case. One can rewrite the baryonic current as 



1 



Jb 



24tt 2 

ie 



32tt 2 

This is manifestly gauge invariant. 



6^ pa Tr[R u R p R a } 
€ ^p* F Tr[<r*(L a +R a )]. (38) 



The first term in the current jg gives a topological 
number associated with tts(SU (2) l-r)- Indeed, the inte- 
gration of it over the space gives the topological winding 
number 



N B 



1**2^ 



Ajk 



Tr[RiRjRk\. 



(39) 



The second term in the current j^ nm is a manifestation 
of the chiral anomaly. Actually, plugging Eqs. (fTTj) and 
([T2]) into Eq. (|35|h one finds the famous tt° — )• 27 term 
by the anomaly 



N c e 2 
~48ir 2 F n 



e^A p F„ p d a ir° + 0{F~ 3 ) 



N c e 2 
12tt 2 F, 



ir E-B + O(F- 3 ), 



(40) 



where the equality holds up to a total derivative. 

Let us next obtain the electric charge coupled to a pho- 
ton fluctuation a M under a background electromagnetic 
field A p. To this end, we expand the gauge field as 



— A^ 



(41) 



Then the WZW action linear in a M gives 

^WZw[fl/x] 



7 ApdviapPv) 



d 4 xe^ upa 



48tt : 



;Tl[R u R p R 



16tt 2 



Wtt 



2Afi(d u Cip)P (T 



(d v A p )P a \ dp 



(42) 



From this, we can read the electromagnetic current 
wzw/H as 



•7em,WZW — ^ 



J, 



em,WZW 



Janm 



-^-[TrlRuRpR^+SedaiApP^ , (43) 



16tt 2 



(d v A p )P a 



(44) 



The current j£ m wzw is gauge- invariant Q • The total 
electromagnetic current is 



Jem — c Jv ' -'em, WZW 



(45) 



This is gauge- invariant, and always conserved at on- shell, 
from the gauge invar iance. 

The total electric charge is given by the spatial integral 
[2lj of the zero-th component of this current ([45]) , 



I 



d 3 x BiP. 



(46) 



Here B\ stands for a background magnetic field and we 
have used e° 2j/e = —e^ k . This is the expression of gauge 
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invariant and conserved electric charge which includes 
an extra term to the well-known Gell-Mann-Nishijima 
formula. The last term is a contribution from the chiral 
anomaly. 

Note that, as we will see below shortly, the surface 
term in Eq. ([43]) does not contribute to the net electric 
charge if the pion mass is non-zero. So, we focus on the 
new last term of Eq. ([46]) coming from Eq. ([44]) . in the 
following of this paper. 



B. The Skyrmion: A nucleon as a topological 
soliton 



Let us find a solution of the classical equations of mo- 
tion derived previously, 



(47) 



We solve these by dealing with the electromagnetic inter- 
action as a perturbation. Then we expand the chiral field 
with respect to the electromagnetic coupling constant e 
as 



Z7 = exp (zr-(/o + A + ---; 

^ = 4 0) + 4 1) + ---- 



(48) 



At the zeroth order the chiral fields and electromag- 
netic fields are decoupled, so that the equations of mo- 
tion are those of the Skyrme model without the electro- 
magnetic interaction and the Maxwell equation without 
a source 

^1° V = -^|^Tr [U - C/t] , dtlF W^ = o, (49) 



with 



(0)/i 



* f F 2 R (0), _ 1 [jR (0) 5 [^(0)^(0)^ (50) 



Jl 



The second equation in Eq. ([49]) is solved by considering 
a constant background magnetic field, say along the xi- 
axis 



1 



ijk F jk = B l 



(51) 



with B l being a constant. 

In order to solve the first equation in Eq. ([49]) . it is 
useful to introduce a dimensionless coordinate 



1 



X U1 



(52) 



In terms of this new coordinate, the Skyrme equation is 
written as 

d M - [RW, [rW",R(°»]]) = ^Tr [C/t _ f/](53) 



with a dimensionless mass in unit of l/(e s F 7r ) 



(54) 



Here and after, we will use this notation. 

Let us make a standard hedgehog (radial) ansatz, for a 
static and topologically non-trivial solution with Nb = 1, 

U (x) = exp (ifo ■ rj = exp (if(r)r-x) , (55) 

Xi = — . (56) 
r 

One can express U in a different fashion as 

U = l 2 cos|/o| + ^777 -rsin|/o| 
\fo\ 

= I2 cos / + ix • r sin /. 
Then we obtain for a static configuration 



(57) 



R 



(0) 

— IJ XiX ■ T 



(58) 



(^ji — x • rxij (s'm2f + 2ix • r sin 2 f^j . 



Putting this into Eq. ([53]) . we obtain the ordinary differ- 
ential equation 



1 2 sin 2 / 



r +¥r f 



sin 2/ fj, 2 _ 1 _ sin / 



m 



^sin/. (59) 



The solution with a unit winding number corresponds to 
lim f(r) = 7T, lim f(r) = 0. (60) 

r— 7-00 r— >-0 

Numerical solutions with different m^'s are given in 
Fig. [2j We adopt the physical pion mass mf ys = 0.263 
which was determined from the mass splitting between 
nucleon and A [l2| . We also show the profile functions for 
= 2.63 and 0.0263 in order to demonstrate a charac- 
teristic property of the profile function. We see that the 
larger (smaller) gives the thinner (fatter) Skyrmion. 
The profile function with = 0.0263 almost coincides 
to that with — 0, see Fig. O 

Asymptotic behavior of / can be found by solving the 
linearized equations of motion for large r 



f"+-f'-%-mlf = 0. 



This is solved by 



/ 



C A 



(61) 



(62) 



where C and A(= Cm^) are constant. For massless pion, 
we find 



/ 



c 



(63) 



6 



rrijj; - 2.63 




m n - 0.263 



nu = 0.0263 



FIG. 2: Profile functions for the hedgehog solution. 



with C ~ 8.638. For massive pion, the asymptotic form 
decays exponentially, 



/ - - e 

r 



(64) 



where A is a constant which depends on m^. For exam- 
ple, we find A ~ 2 for = 0.263. 



III. ANOMALY-INDUCED CHARGE 

We substitute the Skyrme solution to the electromag- 
netic current calculated from the WZW term, to evalu- 
ate the anomaly-induced electric charge of the baryons. 
We here use classical Skyrmions for an illustration first, 
then we move onto quantized Skyrmions, to obtain a for- 
mula for the anomaly-induced charge for baryon quan- 
tum states. The anomaly-induced charge is the last term 
of Eq. (|46|) . Since we are considering a constant mag- 
netic field background, it is enough to see Pi defined in 
Eq. dHTJ). 



A. Classical evaluation 

To evaluate P^, let us first write down the left- and 
right- invariant one-forms as 

Ri = i(r - x^j f'xi + fyx) sin 2f 

± (diX - x^j 1 + i (diX x x^j • r sin 2 /, (65) 

Li = i(r • xj fx{ ± ^ (V • dixj sin 2/ 
± [diX • x^j 1 — i [diX x x^j 

Plugging these into Eq. (f3T|h we have 



sin 2 /- (66) 



Pi = -f f XiX 3 - -(dix 3 )sm2f. 



(67) 





FIG. 3: Nb — 1 Skyrmion solution for = 0. From left to 
right, contours plots of the baryon number density, Pi = ±0.2 
and P3 — ±0.2, respectively. 



The topological charge density, P\ and P3 are shown in 
Fig. [3] at rriir = 0. The induced electric charge densities 
with nonzero (see Fig. 0]) are quite similar to those 
for the massless case in Fig. [3] A tiny difference comes 
form the similar but a little different behaviors in profile 
functions f(r) as shown in Fig. [2j 

The spatial integrations of P\ and P3 become 



with 



d 3 x Pi = 0, / d 3 x P 3 



c 



47T 



r c , (68) 



3(e s P^) 2 

J dr (r 2 f + r sin 2f) . (69) 



Note that the integration variable r —> r/(e s F 7T ) is the 
dimensionless coordinate, so that Co is a dimensionless 
number. This means that the net induced charge is zero 
for B oc (1,0,0), and (0,1,0) whereas it is non-zero for 
B oc (0,0,1) where non-zero correction appears to the 
Gell-Mann-Nishijima formula. The numerical coefficient 
Co can be rewritten as 



c 



[r 2 f]o+ J d r (-2r/±rsin2/). 



(70) 



The first term is surface contribution which becomes non- 
zero only at = due to distinct behavior at large r 
shown in Eq. ([63|) , and the second term expresses a pion- 
cloud effect discussed in Appendix [Bl Computational 
results of Co become 

c =-14.1±C (m^O), (71) 
c = -10.2 (m^ = 0.263), (72) 



where C is asymptotic factor shown in Eq. (|63|) . Com- 
putational results for other values of the pion masses are 
also summarized in Table H 

Finally, we evaluate contributions from the surface 
term in Eq. (|43|) . To this end, we need to compute 



W = 



= -Jd 3 x e^diiAjPk). 



(73) 



Interestingly, plots of the integrand (density) are quite 
similar to those in Fig. [3] if we choose a gauge Ai = 
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FIG. 4: The contour plots of the anomalous charge densities 
of B — 1 Skyrmion for /3 = 0.263 on the cross section by the 
y = plane. The top panel shows —Pi/2 with B\ ^ and 
the bottom panel shows —P3/2 with B3 / 0. The blue lines 
have positive values and the red ones have negative values. 
The black broken lines correspond to zero charge contours. 



£i(0,-2/2,2//2) or A { = B 3 (-y/2,x/2,0). The surface 
term can be evaluated as 



W 



dn 2 



Xie ijk Aj ( f f x k x 3 + 



(d k x 3 )sm2f 



As shown in Eq. (ffTi]h the profile function for 7^ 
is exponentially small at large r, so that this integration 
vanishes. 

Note that in the massless case the asymptotic behavior 
given in Eq. ([62]) leads to W ^ as 



W 



Z2mC 



(74) 



We see that this surface term with the constant C given in 
Eq. ([62]) cannot cancel the last term of Eq. ([46]) . Anyway, 
since the physical pion mass is not zero, we consider the 
new last term in Eq. ([46]) as the anomaly-induced electric 
charge [22j . 



B. Evaluation with quantized Skyrmion 

To evaluate the anomalous current and charge for each 
baryon state, the Skyrmion is quantized as a slowly ro- 
tating soliton. Quantizing the collective coordinates of 
soliton's moduli space G is achieved by the canonical 
quantization of a particle on a manifold G. In the case 
of two flavors, G = SU{2) c± S 3 . We construct the angu- 
lar momentum operators acting on baryon states and the 
harmonic functions on G corresponding to baryon wave 
function. Using these, we evaluate the expectation values 
of the anomalous currents. 

We evaluate the expectation values (jg nm )B and the 
anomalous charge Q^ nm for each baryon state B in the 
presence of a background magnetic field. We show that 
the spatial components of the currents vanish: (jl nm ) b — 
0. We also obtain the anomalous charge Qanm, which is 
given by integrating 0'° nm )s. 



1. Angular momentum operators and spherical harmonics 
on S 3 

Let g £ G be a group element of a group manifold G. 
Then there are operators C a and 7Z a acting on g from left 
and right, respectively, and satisfying the commutation 
relations 

[C a ,Cb] = ifabc£c, [Ro,%] = ifabcHc- (75) 

Here the roman indices correspond to those of the tangent 
space of G, and f abc is the structure constant of the Lie 
algebra of G: [T aj T b ] = if abc T c and Ti(T a T b ) = S ab /2. 
The actions of C a and lZ a on g are 



£ a g = -T a g, C a g 1 = g 1 T a , 
K a g = gT a , K a g~ l = -Tag' 1 . 



(76) 



We restrict ourselves to the case G = SU(2), and thus 
fabc — tabc- This is the case of the two-flavor Skyrmion. 

The operators C a and lZ a are precisely the angu- 
lar momentum operators with respect to the isometry 
SU(2)l x SU(2) R of S 3 . We introduce the scalar spheri- 
cal harmonics on S 3 , Yj m7 n, where J is the same magni- 
tude spins of both SU (2) l and SU (2)#, and m and fh are 
the eigenvalues of their third-components, respectively. 
The actions of the operators on the spherical harmonics 
are 



C±Yj rn , 
Tl±Yjm> 



^(JTm)(J± 



m 



l)*J(m±l)mj 

™>)(J ± ™ + i)^Jm(m±l) , 
CzYjmrh = mYjrnrn, K?,Y Jrnr ~ n = fhY Jrn fh, (77) 



where C± = C\ ± 1C2 and 1Z± =7Zi± HZ2. 

It is convenient to introduce a D-function D ab (g) so as 
to see the relation of Eq. ([76]) and Eq. ([77]) . It is defined 
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by the adjoint action of g, 

gT a g^ = T b D ba (g). 

The action of C a and lZ a on D ab (g) becomes 

C c D ab (g) = ie cad D db (g), 
ft c D ab (g) = -ie cdb D ad (g). 



(78) 



(79) 



With a little more algebra, it can be shown that appro- 
priate linear combinations of D ab (g) precisely give the 
harmonic functions Yjrafa- 

Note that C a and lZ a are respectively the isospin and 
the spin operators for baryon states. See, for instance, 
Ref. [15] for a discussion in three-flavor case. The ad- 
joint action of g to the hedgehog Skyrmion gives U (x) = 
gU (x)g^ = U (x rot ). That is, gx • r#+ = x rot -r. The 
transformation of the unit vector x under the spatial ro- 
tation caused by D(g) is 



x r a ot = D ab (g)x b . 



(80) 



It is natural to identify SU (2)r as the baryon spin, where 
g -> gkR with k R e SU(2) R . 



2. Absence of the spatial anomalous current 



where Cj^^j is a Clebsch-Gordan coefficient of 
SU(2). The wave function of each baryon state is given 
by Yjmrn- Our primary interest is in nucleons (I = J = 
1/2) and A baryons (/ = J = 3/2), but here we can keep 
J arbitrary. We use Eq. ([84]) to evaluate Po projected 
onto each baryon state. 

In Eq. ([83]) . we need to focus only on the last line: 



(85) 



Ow = / l-o£+ + ^l+o£-)weyl 

= Yi- £+ + Yi +0 £_ + V^Yioo. 



It is easily seen that each of the first and the second lines 
in Eq. ([83]) gives no contribution. Integrating Eq. ([85]) by 
using Eq. ([M]) along with baryon states of quantum num- 
bers (J, I 3l 63) for (iso)spin and the third components, we 
obtain 



(O w )i 



dQ 3 
2^ 



J-h)(J + h + i)c£\ 



J(/ 3 +i) 



^ C 10JI 3 



(86) 



To evaluate (jl nm )B in the presence of a background 
magnetic field, we need to focus only on Po owing to the 
index structure of the WZW term. We first write down 
Po and then quantize it. Substituting a slowly rotating 
Skyrmion U(x,t) = g(t)U (x)g^(t) for P (371), we obtain 

P = 2sm(2f)e ac3 D ah (g)x h Tr[r c gg^. (81) 

In the procedure of the canonical quantization, the time- 
derivative part is replaced with the angular momentum 
operator as follows jl6| : 



C a = iATi[r a gg% 
A = ^ / drr 2 sin 2 / 



1 + 4 / 



sin / 



(82) 



Hence Po can be written in terms of raising- and lowering- 
operators, 

Po = sin(2/)[5+(y 1 __r+ + 

-i_(F 1 _ + £++F 1++ £_) 

+ £ 3 (Yi_o£+ + y"i + o£_)]weyi, 

(83) 

where the indices ± in Yjmfh mean ±1, and x± = x 1 ± 
1x2- The Weyl ordering for the operators is understood. 

Integrating a product of three spherical harmonics over 
S 3 gives 



j 2tT 2 ^ Jirni '^ 11 ^ ^ r j2rn 2 rh 2 Yj 3 



/(2J 2 + 1)(2J 3 + 1) 
2Ji + 1 



J 2 rn 2 J 3 m, 3 J 2 fh 2 J 3 rh 3 '> V° / 



where the minus sign appearing in front of S3 is due 
to Eq. (|75j) . This exactly vanishes once values of the 
Clebsch-Gordan coefficients are substituted [17j: 



c 



J 1 

1±1J(7T1) 



+ 1 



' (J± 7 )(J T7 + 1) 
2J(J+1) 



(87) 



Thus we see 0lnm>£ = °- 



3. Anomalous electric charge in baryons 



In j^ nm , rotation of the Skyrmion is encoded in x r 3 ot . 
Hence, it is sufficient to evaluate (x r 3 ot )B, which directly 
leads us to (j® nm )B- Since this part does not contain 
derivatives in time, we simply integrate Eq. ([80]) by using 
Eq. dBU). The result is 



I3S3 



J(J + i) 



x 3 . 



(88) 



Below we will mainly focus on the case J = 1/2 for sim- 
plicity. However, thanks to Eq. ([88]) , it is straightfor- 
ward to consider higher-spin cases. For instance, this 
gives (x r 3 ot ) N = -4/3S3X3/3 for nucleons, and (^ ot ) A = 
— 41363X3/15 for A baryons. 

Let us calculate the total electric charge from the 
anomalous effect. The matrix elements of P, are eval- 
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uated by applying Eq. 

(P ): 



)ls 3 = o, 



,p \N a=l,2 
\^a/I 3 ,S 3 



hs 3 f 



(89) 

sin(2/)\ _ 

l )x a x 3 , (90) 



hs 3 



r 



2r 
sin(2/) 
2r 



sin(2/) 



2r 



(91) 



The anomalous charge density under a constant magnetic 
field B is indeed induced in nucleons: 



■0 \N _ ie N c p /pvJV 

tl3,S 3 - 487T 2 ^^3,53" 



wanm/ i 



(92) 
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The integration of (^i)^ s 3 over the whole space yields 



167TZ 



(4I 3 S 3 )co (i 



1,2), 
3), 



where Co = J dr{r 2 f + r sin(2/)}. Numerical value of Co 
will be shown in Table [I] for several pion masses. From 
Eq. ([92]) , we obtain the anomalous charge for nucleons 



✓ anm 



4e 2 N c 
27tt 



CqB 3 



(93) 



In this final expression we restored the rescaling factor 
e s F n by a dimensional counting. Equation ([93]) shows 
that an electric charge is actually induced by the anoma- 
lous effect even for a neutron. 

As seen from Eq. ([88]) . dividing the result in Eq. ([93]) 
by a factor of 5 gives the anomalous charge of A baryons. 

The plot of the charge density of (j® nm ) for the quan- 
tized Skyrmion shows exactly the same as Fig. [3] For 
a magnetic field along x 1 direction, the charge density 
plot is symmetric, thus the total charge vanishes. How- 
ever, obviously multipoles, in particular a quadrupole, 
may show up. In the next section, we calculate multi- 
poles in (j° nm ). 



IV. MULTIPOLE MOMENTS OF ANOMALOUS 
CHARGES AND PION-MASS DEPENDENCE 

When one regards charged baryons as point-like parti- 
cles, multipole moments are suitable physical quantities 
to describe an original charge distribution. In this sec- 
tion, we extend the calculations to the higher multipole 
moments due to the anomalous- charge distributions, and 
estimate pion-mass dependence of the anomalous charges 
and the multipole moments. 

First, we can easily find that the dipole moment due 
to the anomalous charge vanishes: 



Di 



d 3 x Xi (f anm ) N = = 1,2,3). 



(94) 



FIG. 5: Behavior of the Skyrmion profile functions, /(r), for 
several pion masses. 



10 z 



10 1 



10 u 



10" 



Fit: -6.0(7) /4 97(5) 



+ + 



+ \ 
+ 



-c + 



0.1 



10 



100 




100 



FIG. 6: Results of the numerical coefficients, — Co (top) and 
— C2 (bottom), as a function of /ra£ hys in log- log scale. The 
dotted lines mean fit results by a function A/m™ with A and 
n being free parameters in the fit ranges shown by the arrows 
on the figures. Numerical results of the fit are also shown on 
the figures. 
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On the other hand, the quadrupole moment: 

Qij = J d 3 x (3x iXj - r 2 5 tj ) {j° an J N , (95) 
is calculated for the nucleon as 

= e l^(hS?)Qijk , eB ^ A c 2, (96) 

(-2S k 3 3£fcA 
-25 ks 3S k2 , (97) 

where the numerical coefficient, 

c 2 = J dr[2r 4 /'-r 3 sin(2/)] , (98) 

is shown in Tab. [I] for several pion masses. This means 
that the leading multipole due to the anomalous con- 
tribution is the quadrupole moment. We note that the 
quadrupole is induced in response to all directions of the 
external magnetic fields, although the anomalous charge 
is induced only by B3 (see Eq. ([93]) ). 

In order to extract the pion-mass dependence of the 
anomalous charge and the quadrupole moment, we cal- 
culate the Skyrmion profile function f(r) for wide pion- 
mass range (0.1 < m 7r /mP hys < 100). Behavior of f(r) 
for several pion masses is shown in Fig. [5j where the solid 
line is f(r) at = mf ys . We find that the wave func- 
tion shrinks with the pion mass increasing. 

Since the pion mass dependence of the anomalous 
charge and the quadrupole moment appears in the nu- 
merical coefficients, Co and C2, via r integration with 
/(rjraTr), we focus on these coefficients. Figure [6] shows 
results of — Co (top) and — C2 (bottom) as a function of 
m^/m^ ys in log- log scale. Numerical values of Co and C2 
are also summarized in Tab. HI We can see that Co be- 
comes almost plateau at small pion-mass (m^ / m^ hys < 
1), whereas that decreases linearly at large pion-mass 
(ra^/raP hys > 10). We fit the results by a function, 
A/m™, with A and n being free parameters, and obtain 

c 6.0(7)/ra°' 97(5) in a range of 30 < m./mfy 8 < 

100, shown by the dashed line on Fig. [6] (top). This im- 
plies that the pion-mass dependence of the anomalous 
charge is Q a nm oc at large pion-mass. 

On the other hand, C2 behaves linearly for all pion-mass 
region. We also fit the results by A/m™, and obtain C2 ~ 

-127(15)/ra^ 06(7) in a range of 20 < m^/m^ 8 < 100, 
shown by the dashed line on Fig. [6] (bottom). Although 
the fit is performed at large pion-mass region, almost 
all results of C2 is located around the dashed line. This 
implies that the quadrupole moment due to the anomaly 
behaves as Qij oc 1/m^.. 

Note that we have evaluated j^ nm which is only a part 
of the total electromagnetic current. (23| 





Co 


C2 





-1.41X10 1 


— 00 


0.1 


-1.37X10 1 


-4.90xl0 5 


0.125 


-1.37X10 1 


-2.99xl0 5 


0.25 


-1.34X10 1 


-6.68xl0 4 


0.5 


-1.23X10 1 


-1.50xl0 4 


1 


-1.02X10 1 


-3.27xl0 3 


2 


-7.32 


-6.91xl0 2 


3 


-5.67 


-2.76xl0 2 


4 


-4.62 


-1.44xl0 2 


5 


-3.88 


-8.63X10 1 


1 n 


— 9 90 


— 1 £Q v 1 D 1 


20 


-1.15 


-4.14 


30 


-7.97xl0 _1 


-1.74 


50 


-4.99xl0 _1 


-6.90x10" 


100 


-2.38xl0 _1 


-1.35x10" 



TABLE I: Numerical results of the coefficients Co and ci. We 
neglect the surface contribution at = shown in Eq. ([7T|) 
due to distinct behavior of the Skyrmion profile function f(r) 
at r —> 00, i.e. f(r) ~ 1/t 2 , discussed in Sec. Ill Bl 

V. NO CONTRIBUTION FROM OTHER 
CORRECTIONS 

In this section, we study other effects of the back- 
ground magnetic field to the the electric charge of the 
nucleon. Our aim is to show that the anomaly-induced 
electric charge is not cancelled by the other electromag- 
netic effects. 

The total electric charge is written as a modified Gell- 
Mann-Nishijima formula, 

Q e = e/ 3 + ^ + ^p. (99) 

The first term stems from the electromagnetic current 
in the original Skyrme model. The second term is due 
to the baryon number coupling to the electromagnetic 
potential. The last one is the anomaly-induced electric 
charge which is nonzero only when we have the back- 
ground electromagnetic field. 

We are working with the perturbative expansion with 
respect to the background magnetic field eB. What we 
found for the anomaly-induced charge is 

Qanm = D{e 2 B). (100) 

In the total charge formula ([99]) . the second term is due 
to the baryon charge which is a topological charge for 
the Skyrme model, thus not corrected by the background 
magnetic field. On the other hand, the first term can be 
corrected in the presence of B. If a correction of the order 
0(eB) appears from the ^3 term in the charge formula, 
then it may possibly cancel our anomaly-induced charge. 
In the following, we shall present an argument showing 
that there is no such correction of 0(eB) to the ^3 term. 
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First, let us examine if there is a correction to the 
electromagnetic current itself in the Skyrme model. One 
may naively think that, since the Skyrme solution it- 
self is corrected by the background electromagnetic field, 
the current may also be corrected. However, this is not 
the case for the Skyrmion. The reason is that for the 
Skyrmion, the electromagnetic ?7(1) is identical to a part 
of the isospin, and the action itself has the isospin struc- 
ture from the first place. In fact, the relevant I3 is indeed 
expressed by a part of the flavor SU (2) rotation, and thus 
the current Is is universally expressed as 



a 



_d_ 

da 3 



as 



d 
da 



ai 



_d_ 

da 2 



CL2 



_d_ 

dai 



(101) 



Here, there is no room for eB to show up, thus we can 
safely use this expression for I 3 in the electric charge 
formula ([99]) . 

Then, the issue is whether the expectation value of Is 
in the background magnetic field is corrected or not. The 
background magnetic field modifies the wave function of 
the Skyrmion, so in principle this corrected wave function 
may give a correction to (Is), which is of importance for 
us. We shall show in the following that there is no such 
correction at O(eB). 

To proceed, we need to know how the Skyrme moduli 
wave function is corrected. Due to the background mag- 
netic field, there appears a potential in the moduli space, 
then two of the moduli parameters are lifted to become 
pseudo- moduli parameters. The corrected quantum me- 
chanics of the moduli and the pseudo-moduli is written 



3 



(102) 



where the potential of the quantum mechanics is of the 
form V(a) = ((ai) 2 + (a 2 ) 2 ) V(a). The function V(a) is 
a polynomial of (with a finite order), with just numer- 
ical coefficients. The potential V(a) breaks the SU(2) 
symmetry of the system down to the diagonal U(l). 

We briefly explain how to derive the form ([102)) of the 
induced potential. In the Skyrme model, the electromag- 
netic interaction enters as 

R„ = D^UU\ D^U = d^ ieA^q, U}. (103) 

For a background magnetic field Bs, we consider A\ = 
—Bsx 2 , thus among the electromagnetic contribution 
appears only in R\ as 

£1 =d l UU 1[ + SR U (104) 

SRx = -ieBsx 2 [q,GU G^]GU^Gl (105) 

Here we have already substituted the Skyrme solution. 
Assuming that G is dependent on t and plugging this 
R into the action, we obtain a correction to the Skyrme 
Lagrangian at O(eB) as 

SL = ^TrlRtfRi] + ^Tr[i? M , R^R?, (106) 



Since we know that the Skyrme solution has the partic- 
ular dimension dependence x — )• x/(e s F 7r ), we obtain 



SS= fd 4 xSL= fdt^-eBsV(a), 

J J e s-^7T 



(107) 



where V(a) is a polynomial in a^, with only dimensionless 
numerical coefficients. This is nothing but the potential 
in Eq. (|102p . As the potential V(a) should vanish when 
G corresponds to the electromagnetic direction, i.e. the 
rs direction, V(a) is proportional to (ai) 2 + (a 2 ) 2 . 

With the potential, the Skyrme wave function ij)(a) is 
modified. We may apply a well-known perturbation tech- 
nique for quantum mechanics, and obtain the corrected 
nucleon wave function as 



1/2) = \l = l/2)o 



3*3 £ 



Vl=l/2, l=n+\/2 | 
" E l=1/2 ~ E[ =n+1/2 



0{{eBf). 



n+l/2) 
(108) 



Here ^=1/2, i=n+i/2 1S the matrix element of the oper- 
ator V appearing in the quantum mechanics (jlQ2j) . and 
the state with subscript is the one without the pertur- 
bation. In the current case the states have degenerate 
energy, but the expression above is universal. 

Now, using this corrected wave function, we evaluate 
the expectation value of Is- Since we have 



o(Z = n+l/2|/ 3 |Z = l/2) = 
for n > 1, we obtain 

(I 3 ) = (I 3 ) + O((eB)% 



(109) 



(110) 



where (13)0 is the third component of the isospin of 
the leading (uncorrected) order wave function. There- 
fore, the electromagnetic correction to the charge formula 
starts at (D((eB) 2 ), which is at higher order compared to 
the anomaly-induced charge Q a nm- This means that our 
anomaly-induced charge Qanm is the leading-order cor- 
rection of 0{eB\ and cannot be cancelled by the other 
effect of the background magnetic field. 

Here we have presented an argument that the total 
induced charge due to the anomaly is not cancelled by 
other corrections due to the magnetic field. This argu- 
ment may be reenforced and supplemented by an explicit 
computation of a back-reaction to the Skyrme configura- 
tion itself, from the magnetic field. The calculation of 
the back-reaction is quite complicated, so we leave it to 
our future work. 



VI. HIGHER-CHARGE SKYRMIONS 

In this section, we study classical higher-charge 
Skyrmions and the anomaly-induced charges. To this 
end, we will utilize the so-called rational map ansatz [l9| 
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which is a reasonable method giving a good approxima- 
tion. In this section we use another notation based on 
a standard textbook [20]. A main difference from the 
previous sections is the dimensionless coordinate 



M e F 



(111) 



Let us first give a brief review on the rational map 
ansatz. A solution of the Skyrmion U(x) with U(x) — )• 1 
as \x\ — » oo gives a map from R 3 + {oo} ~ S 3 to 
SU(2) ~ S* 3 . The map is characterized by the homo- 
topy group 7Ts(SU(2)) = Z. More explicitly, it can be 
expressed as 



U(x) = exp (if B (r)r • n(#, < 



(112) 



where {f Bl n} is a coordinate of SU (2) under a constraint 
Ib £ [0,tt] an d \n\ = 1: Namely, we decompose SU(2) 
into I[o,tt] x S 2 . The parameters (r, 0, 0) are standard 
spherical coordinates on the space R 3 c± R>o x 5 2 . In 
order to get the map of degree E Z, we assume that 
f B is a one-to-one map from R>o — > /[o j7 r]- Then n(#, 0) 
should give a map 5 2 — >> 5 2 of degree N B . Let us intro- 
duce the stereographic projection which is useful to find 
the generic map of degree N B , 



By making use of the following pull-back 



i + N 2 x: 

1 + \w\ 2 



dw 



(119) 



it is easy to change the integral area over S 2 to the target 
space of the rational map S 2 as 



r.h.s of Eq.Gni! = - J ^ sin2 ^ 



drdSlo, 



2N f 



sin f B df B = N B , 



(120) 



where we have used J dVt w = 4ttN b . The Skyrmion en- 
ergy in the F n /(4e s ) energy unit and 2/(e s F 7T ) length 
unit with the rational map ansatz is given by 



E = 4tt / dr 
Jo 



+J 



where we have introduced 



r 2 f B 2 +2N B (f B 2 + 1) sin 2 f B 
sin f B 



+ 8m 2 T (l-cos/ B ) 



. (121) 



z(6,(f>) = e # tan-. 



(113) 



For example, the N B = 1 hedgehog ansatz (the one-to- 
one map from S 2 to S 2 ) can be expressed by 



n(z, z) 



z 1 



1 



1 



1 



(114) 



This can be easily extended to a iVe-to-one map from S* 2 
to S 2 by replacing z with any rational maps w(z) : C to 



1 47T 



1 + M 



dw 



dn z . 



(122) 



In order to find the best approximation, we need to 
seek an appropriate rational map w(z). We should choose 
w in such a way that / is minimized. Though this is 
not easy task, by using a numerical method, the rational 
maps w(z) for several N B were found in Ref. [19]. For 
instance, the following rational maps for N B 
are known as 



1,2,--. 



w(z) 



<?(*)' 



(115) 



Here P(z) and Q(z) are holomorphic functions in z and 
we set N B = maxjdegP, deg Q}. Thus, we obtain the 
map from S 2 to S 2 of degree N B , 



w + w . w — w 1 — \w\ z 



1 



w 



1 



w ■ 



1 



w 



(116) 



Plugging this into Eq. (|112jl . we reach the map from S 3 
to SU(2) of degree N B . This is called the rational map 
ansatz. 

The baryon number can be expressed by 



=-[Ik. 

' J 2** 



f B /sin/fl l + |z| 



dw 



r 2 drdn z , (117) 



r 1 + \w\ 

where dft z is the usual area element on a 2-sphere 

2idzdz 



dtt z 



(1 + M 2 ) 2 



sin OdOdr. 



(118) 



Wi 

w 2 
w 3 

W4 

W 5 
Wq 
Wj 

w 8 



y/3iz 



VSiz 2 - 1 ' 
z 4 + 2V3iz 2 + 1 
z 4 - 2\[liz 2 + 1 ' 
z (z 4 + bz 2 + a) 



6z 2 
- zc 



z 2 (icz 4 + l)' 
7 - 7z 5 - 7z 2 - 1 



z 7 + 7z 5 - 7z 2 ■ 
z 2 (dz 6 + l)' 



1 



(123) 
(124) 

(125) 
(126) 

(127) 
(128) 
(129) 
(130) 



with a 3.07, b 3.94, c 0.16 and d 0.14. 

The last task is to determine the profile function f B . 
Because no analytic solutions have been known, we need 
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r 



5 10 15 20 



FIG. 7: The profile functions f B =2 for N B = 2. 



to solve equations of motion numerically, 

(i + ^*ff B y B+ 2 -f B 

N B sm2f B ( fl2 I sm 2 f B \ 

+ r 2 [Jb 1 Nb r2 J 

-4m2 s in/ B = 0, (131) 

with the boundary condition /#(()) = tt and /s(oo) = 0. 
We show several numerical solutions for Nb = 2 with 
different pion masses in Fig. 

Now we are ready to evaluate the anomaly induced 
electric charge from Eqs. ([37]) and ([46]) . As before, what 
we need is only Pi which can be obtained from 

Ri = i(r -n)f' B Xi + -(t- din) sin2f B 

+ {{dS ■ n)l + i(difi x ft) • r] sin 2 f B , (132) 

L { = i(r ■ n)f' B Xi + -(r • ^n) sin2/ B 

+ {(ftn • n)l — z(^n x n) • r} sin 2 (133) 
By plugging these into Eq. ([37|) , we get 

^ = -/>3^ - ^(a,n 3 )sin2/ B . (134) 

Note that, as expected, replacement 77,3 with X3 gives us 
N B = 1 hedgehog solution. 

The induced charge densities for the N B = 2 solu- 
tion are shown in Fig. [51 As one can see, the N B = 1 
and N B = 2 charge distributions are quite similar, even 
though the baryon charge distributions are totally differ- 
ent. However, one can find differences if paying attention 
to the detail structures. As can be seen in Figs. [4] and 
[9j the N B = 2 densities are fatter than those of N B = 1. 
Also the N B = 2 configuration has an internal structure. 

The anomaly induced charges of the classical 
Skyrmions with N B = 1, 2, • • • ,8 under a constant back- 




FIG. 8: Nb = 2 Skyrmion solution. From left to right, the 
baryon number density, Pi and P3, respectively. 




-^-20 2 4 



FIG. 9: The contour plots of the anomalous charge densities of 
B = 2 Skyrmion on the cross section by the y = plane. The 
top panel shows — Pi/ 2 with B\ 7^ and the bottom panel 
shows — P3/2 with B3 / 0. The blue lines have positive values 
and the red ones have negative values. The black broken lines 
correspond to zero charge contours. 



ground magnetic field B = (0, 0, B3) 

e 2 

°— ^ = l^(e s F^ B ^ (135) 

c = 4 J d 3 x P 3 , (136) 

are summarized in the Table [III The pre-factor 4 is 



N B 


1 


2 


3 4 


5 6 7 8 


CO 


-43.2 


-105 


-60.3 0.00 


-13.3 28.7 0.00 -11.6 



TABLE II: The anomaly induced charge of N B = 1, 2, • • • ,8 
Skyrmions under a constant magnetic field background. The 
dimensionless pion mass is chosen to be — 0.263. 



needed because of the dimensionless coordinate —> 
2x fI /e s F 7r . Note that Co and Co for Nb = 1 are related 
by Co = 47rco/3. We find that the classical anomaly- 
induced charge is not proportional to the baryon charge 
Nb. It is intriguing that Nb = 4 and 7 Skyrmions have 
zero induced electric charge. From the values given in 
the Table [III we observe that higher-charge Skyrmions 
tend to cancel the total induced charge. A natural rea- 
son for this cancellation is as follows. Each Skyrmion 
has a classical orientation in spin and isospin space, and 
to form a bound state of the Skyrmions the orientations 
should be arranged to cancel each other. Our formula 
of the anomaly-induced charge depends on the signs of 
the quantum spin and isospin, so, accordingly the to- 
tal anomaly-induced charge would tend to cancel each 
other. Although we have not performed quantization of 
the higher-charge Skyrmions, we expect that this cancel- 
lation should occur even at the quantized level. 

Let us finally display the baryon number densities 
and anomaly induced electric charges of the Skyrmions 
with N B = 3,4,- ,8 and N B = 17 [10], see Fig. [TOJ 
The anomaly-induced charge densities exhibit amusing 
shapes. Possible interpretation of the shape is an open 
question. 



VII. CONCLUSION AND DISCUSSION 

We have evaluated the gauged WZW term for quan- 
tized Skyrmions under a background magnetic field in 
the expansion of the electromagnetic coupling constant. 
We have found that there is an anomaly-induced charge 
structure due to the gauged WZW term. The detailed 
analysis of the total induced charge suggests that the 
pion cloud of the baryons can induce a net charge. The 
magnitude of the induced charge structure is G(e 2 B), so 
it is quite small except for the case with strong magnetic 
field background. 

We have calculated the anomaly-induced electric 
charge for any baryons which appear as quantum excita- 
tions of Skyrmions (Section IIII)) . The induced charge is 
non-vanishing when the magnetic field is present along 
the axis of the quantization of the spin and the isospin 
of the baryon. In Section [Vj we argued that this in- 
duced charge may not be cancelled by other possible 
electromagnetic corrections to the Skyrmion, although 
a complete verification may need an explicit calculation 
of the back-reaction of the Skyrmion solution in the mag- 
netic field. We further examined explicitly the anomaly- 
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FIG. 10: Higher-charge Skyrmion solution. From left to right, 
the baryon number density, Pi and P3, respectively. Nb = 
3, 4, • • • ,8 and 17 are shown from top to bottom. 

induced quadrupole moment (Section HVj) and also the 
cases with multi-baryons (Section IVIj) . 

It is nontrivial that an additional electric charge of 
baryons is generated in magnetic fields. It may have 
an observable effect on physics related neutron stars and 
heavy ion collisions [lOj. 
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Finally we discuss possible origin of the anomaly- 
induced charge. One may wonder if the constant mag- 
netic field may be too artificial and it might be a reason 
for the anomaly-induced charge. In Appendix [Aj we con- 
sidered a magnetic field generated by a circular electric 
current, and we found that the calculated induced charge 
is again nonzero. It suggests that the induced charge is 
not an artifact of the everywhere- const ant magnetic field. 

Then what is the origin of the additional charge? A 
good indication comes from the peculiar property of 
baryons. As shown in Appendix [Bj we found that the 
total induced charge is due to the multi-pion effect in 
the nonlinear sigma model. As the Skyrmion profile ex- 
tends to the spatial infinity, the charge distribution also 
has a tail which elongates to the spatial infinity. This 
would be the origin of the generation of the additional 
electric charge. Obviously, if quarks are completely con- 
fined, the total charge of any baryon should be quantized 
to be a half-integer. However, in reality, any baryon is 
surrounded by a pion cloud, which means that quark- 
antiquark pair can percolate out of the mean volume 
of the baryon. We can interpret our anomaly-induced 
charge as an effective charge carried by the pion cloud 
surrounding the baryon. To make sure our interpreta- 
tion, it is important to calculate the complete effect of 
the magnetic field, i.e. the back-reaction to the Skyrmion 
profile due to the magnetic field. 

The anomaly-induced charge may appear to violate 
the charge conservation. In general, any electromagnetic 
current should be conserved at on-shell when the total 
system is gauge- invariant, and this applies surely to our 
case. However, we considered in our paper only a static 
situation, so we have not considered the situation where 
one turns on the magnetic field gradually from zero to a 
nonzero value, in a time-dependent manner. To under- 
stand the origin of the additional charge concretely, one 
needs to calculate the back-reaction and also the time- 
dependent magnetic fields. We leave it to our future 
work. 
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Appendix A: Anomaly induced charge in circular 
electric-current 



In the above argument of the anomaly induced charge, 
we have assumed a uniform external magnetic-field. 
However the magnetic field should be always closed un- 
less the magnetic monopole appears. In this appendix, 
we consider the anomaly induced charge in the external 
magnetic field generated by a circular electric-current, 
which is instructive for us because the magnetic field is 
closed with finite circular radius, whereas that becomes 
uniform when a radius of the circular electric-current be- 
comes infinity. Here we suppose that an electric field is 
not induced by the electric current. We will show that 
the anomalous charge is induced in the circular electric- 
current even with the finite radius. 

Let us suppose the circular electric-current density 
with a radius a on x?/-plane as, 



j(r) = J ^5(z)5(y/rf+y* - <0(- sin C, cos C, 0), (Al) 

Z7T 

where we assume that magnitude of the electric-current 
is proportional to the radius a. The magnetic field gen- 
erated by the electric-current density can be given by, 



B(r) = f rotfdr^L 
4tt ./ r - r 



(A2) 



with \i being a magnetic permeability. For simplicity, we 
omit the factor fi/Air in the following. One can easily see 
that, in the large radius limit (a — >• oo), the magnetic 
field becomes 



£i=5 2 =0, B 3 =j . 



(A3) 



This is the same situation with the uniform external 
magnetic-field to z-direction. 

When the nucleon is located at the center of the circu- 
lar electric-current, the anomalous charge in the external 
magnetic field is given by an integration of (j anm )jj s , 
shown in Eq. (|92|h over the whole space, 



Q^= i -^fd 3 xB i (P i )l S3 . (A4) 

Notice that the magnetic field is also functions of the 
coordinate variables. Performing the integration over the 
whole angular-space, we can separate three components 
of the anomalous charge: 



xy(r) = J dVt 2 xiXzBi = J d 

aM=/ 



dfl 2 B 3 , p Zf2 (r) = / dn 2 x 2 3 B : 



where Bi = Bi/jo. Then the anomalous charge can be 
rewritten as 

_ 4eN r cjq 

\^ b V 71 tp \ 2 V c *y + c zA + c z,2) , 



27tt 
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FIG. 11: Results of p xy , p Zj ± and p 2 ,2 as a function of r in the 
case of a = 1. 
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FIG. 12: Magnitude of c Xi , + c z ,i + c z ,2 as a function of the 
radius a of the circular electric-current for ra^ /ra£ hys = 0.5, 
1.0 and 2.0. The arrows on the right side denote magnitude 
of c . 



with the numerical coefficients: 

Cx y = ^ / dr \- 2r2 f' ~ rsin ( 2 /)] Pxy(r), 
3 f°° 

c z ,i = — / drrsm(2f)p Zil (r), 
3 f°° 

c z , 2 = ^J dr i 2r2 f ~ rsm(2f)] Pz*{r). 

Namely, we denote c xy (c 2 ,i and c z ^) as component (s) 
of the anomalous charge induced by B x and B y (B z ) 
generated by the circular electric-current. With these 
definitions, one can also show that c xy + c Zi i + c Zi 2 = Co 
at large radius limit (a — >• oo). 

Fig. fTTI shows results of p^, and p^,2 as a function 
of r in the case of a = 1. We find that shows small 
but finite value with a peak at r = a, which implies that 
the anomalous charge is induced by not only B z but also 
B x and B y . p Zj i shows intrinsic behavior: it becomes 



constant at r < a, whereas it vanishes at r > a. p z ^ 
shows smooth behavior without any singularity at r = a. 

Magnitude of the coefficients, c xy + c Zy i + c Zy 2, is shown 
in Fig. [12] as a function of the radius a for m^/m^ ys = 
0.5, 1.0 and 2.0. The arrows on the right side denote mag- 
nitude of Co- The coefficients, c xy + c Zy i + c Zy 2, have finite 
value even with finite radius which implies the anomaly 
induced charge by the closed magnetic field, and converge 
to Co at large radius. It is also found that the coefficients 
shows minimum values at some radius, e.g. a ~ 12 for 
/ m^ hys = 1.0. This may be understood as follows: in 
the case of a ^ rn^ 1 , the Skyrmion feels similar magnetic 
field to the uniform one to z-direction, which induces sim- 
ilar anomalous charge, i.e. c Zj i + c Zy 2 ~ Co- Furthermore, 
since the anomalous charge is also induced by B x and 
By at finite radius discussed above, there is finite con- 
tribution, \c xy \ > 0. This extra contribution gives larger 
induced charge than that induced by the uniform mag- 
netic field. 

We also calculate the multipole moment due to the 
anomaly, and find that the results are similar to the case 
of the uniform magnetic field: the dipole moment van- 
ishes, whereas the quadrupole moment Qij shows finite 
values only for diagonal parts (i = j). 



Appendix B: Multi-pion effect and comparison with 
point-particle picture 



Here we argue that the anomaly-induced electric 
charge is due to the pion cloud which exists around any 
baryon. The pion "cloud," which is the multi-pion effect, 
in the anomaly term is simply the terms with higher pow- 
ers in the tt field. The anomaly term in the gauged WZW 
term Swzw can be expanded as 



#wzw ~ J d A x A0B3P3 

~ J d A x Ti[T 3 U^dU]A B 3 

~ J d A x [diro + ttttOtt + • • • ] A B 3 . (Bl) 



The first term is responsible for the famous ttq —> 27 
interaction, while the remaining terms are the pion cloud. 

In the following, we shall see that, only with the first 
term, the anomaly-induced total charge Q a nm vanishes. 
So, our anomaly-induced total charge is due to the pion 
cloud. 

For the Skyrme solution, we have tto ~ f(r)x 3l so the 
total electric charge induced by the first term in Eq. (|B1|) 
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is proportional to 



d 3 x d 3 ir = / d 3 xd 3 (f(r)x 3 ) 



= 2tt 


^ r 2 sin OdrdO 


"(/'-{) cos 2 * 


+ Z 1 

r . 


47T 

= T 


POO 

/ dr (r 2 /' + 


-2r/) 




47T 

Y 




(B2) 



The last expression vanishes for nonzero pion mass, be- 
cause f(r) decays exponentially at large r, and /(0) is 
finite. So, the anomaly-induced total charge vanishes if 
one use only the single-pion term in the anomaly term 

m - n 

It was discussed in [18| that the anomaly-induced total 

charge of nucleon vanishes, by using a generic argument 
without using the specific Skyrme model. The argument 
[l8| uses only the single-pion term, so our result is con- 
sistent with it. 

Before going to the multi-pion term, we note that, in 
the chiral limit where the pion mass vanishes, the last 
expression is nonzero, since / ~ r -2 at large r (see 
So, in the chiral limit, contribution which comes from the 
single-pion term is nonzero. This is again consistent with 
the discussion in [18[ where the pion momentum is ne- 
glected compared to the pion mass to show the vanishing 
total charge. Note that this discussion on the chiral limit 
is suggestive but not so firm since various observables in 
the Skyrme model diverges in the chiral limit. 

Now, let us evaluate the multi-pion term in Eq. (jBl|) . 
The representative 3-pion term is evaluated as 

J d 3 x TTTrdir ~ J rf(r) 3 dr, (B3) 

which is nonzero for any pion mass. Therefore, we con- 
clude that our anomaly-induced total charge is due to 



the multi-pion effect. 

The point-particle picture of [l8j shows that the 
quadrupole moment is induced as a leading moment. So 
let us compare conclusion of the Skyrmion with that of 
the point-particle picture. 

The anomaly-induced quadrupole moment has been 
written in the point-particle picture as [l8| : 



N c a g A 
6tt (f^ruir) 2 



A/Vr 3 A^', (B4) 



where a = e 2 /47r, and qa and N are the axial cou- 
pling constant and the nucleon wave function, respec- 
tively. In the Skyrmion, the quadrupole moment due to 
the anomaly is given in Eq. ([96]) , where the pion-mass 
dependence of the coefficient becomes, 



A 



C2 ^(m,/e s F n r (B5) 
Using the formula of gA = —7rD/3e 2 in the Skyrme 
model, where D is the numerical coefficients of r inte- 
gration including the pion profile function we can 
rewrite the quadrupole moment with familiar physical 
observables as, 



Qi 



SN r a A 



I3S3 



9A 



457r D (F^m^) 2 



(B6) 



We have checked that the numerical coefficient D does 
not show singular dependence on m n . This implies that 
the pion-mass dependence of the quadrupole moment is 
qualitatively consistent between the point-particle pic- 
ture and the Skyrme picture. 

As a consequence of the comparison, we find no con- 
tradiction between the point-particle picture and the 
Skyrme picture. For further understanding of the 
anomaly-induced charge, calculations of the multi-pion 
effect in the point-particle picture are required. 
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